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Abstract 



A precise determination of \V u b\ can be obtained exploiting the sum rule for 
inclusive charmless semileptonic 5-meson decays. The sum rule is derived on 
the basis of light-cone expansion and 6-flavored quantum number conserva- 
tion. The sum rule does not receive any perturbative QCD correction. In 
this determination of IV^I, there is no perturbative QCD uncertainty, while 
the dominant hadronic uncertainty is avoided. Moreover, this method is not 
only theoretically quite clean, but experimentally also very efficient in the 
discrimination between B — > X u iu signal and B — > X c iu background. The 
sum rule requires measuring the lepton pair spectrum. We analyze the lepton 
pair spectrum, including the leading perturbative and nonperturbative QCD 
corrections. 
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I. INTRODUCTION 



The measurement of the Cabibbo-Kobayashi-Maskawa (CKM) matrix [[I|] element \V u b\ 
is currently one of the important goals of B physics. Standard model predictions employ the 
fundamental parameter \V u b\ as input. As the experiments at the B factories are starting to 
take data, a precise determination of \ V u b \ is increasingly vital for testing the standard model. 
For example, the standard model prediction for the CP-violating asymmetry in B — > J/ipKs 
decays depends on the value of \V u b\. Evidence for the CP- violating asymmetry has already 
been provided by the CDF Collaboration 0. As the dedicated experiments are expected 
to achieve much greater precision, an increase in the accuracy of \V u b\ is highly desirable 
in order to test whether the complex phase of the CKM matrix is the only source of CP 
violation. On the other hand, despite empirical successes, the standard model is definitely 
not the final theory of particle physics. There remain many outstanding issues: electroweak 
symmetry breaking, fermion masses and mixing, CP violation, replication of families, etc. 
A precise determination of \V u b\ is also an important link in pursuit of deeper principles to 
extend the standard model. 

\V u b\ can be determined from both inclusive and exclusive charmless semileptonic B me- 
son decays. The precision of the determination is limited by experimental and theoretical 
difficulties. Charmless semileptonic B decays are a rare process. The main experimental 
difficulty in observing signals from b — > u£u processes is the very large background due to 
b — > civ whose branching fraction is two orders of magnitude larger. On the theoretical 
side, calculations are needed in order to relate the measured quantity to \V u b\ to extract a 
value from the data. Strong interaction effects complicate the calculations, causing theo- 
retical uncertainties. The uncertainties result not just from nonperturbative QCD which is 
notoriously difficult to calculate, but also from perturbative QCD although in priciple it is 
calculable in perturbation theory. With the developments in nonperturbative techniques, 
the uncertainty from perturbative QCD due to the truncation of perturbation series in prac- 
tical calculations can be comparable to that from nonperturbative QCD, indicating that 
these two sources of uncertainty can be equally important. 

Two approaches have been developed for a QCD treatment of inclusive decays of heavy 
hadrons. One called the heavy quark expansion approach || is based on the operator prod- 
uct expansion. One called the light-cone approach |§-[^] is based on the light-cone expansion, 
using the notion and method of deep inelastic lepton-hadron scattering. In both approaches 
the heavy quark effective theory || is exploited as a limit of QCD. There is important 
distinction between the two approaches. The former is a short-distance expansion in local 
operators of increasing dimension, while the latter is a non-local light-cone expansion in 
matrix elements of increasing twist. Moreover, the heavy quark expansion approach invokes 
quark-hadron duality, which assumes rates evaluated at the parton level to be equal to ob- 
servable rates summed over a sufficient number of hadronic channels. Quark-hadron duality 
cannot be exact The heavy quark expansion approach has to use quark kinematics 

and cannot account for the rate due to the extension of phase space from the parton level 
to the hadron level. The rate missing JTTJ] is just a manifestation of duality violation. In 
contrast, the light-cone approach does not rely on the assumption of quark-hadron duality. 
Rates are evaluated in the light-cone approach using physical phase space at the hadron 
level. The light-cone approach was criticized in Refs. fl2fl . However, it has been shown (11 
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that that criticism is mistaken. The validity of the light-cone approach has been tested || 
experimentally on the lepton energy spectrum in inclusive semileptonic decays of B mesons. 
It is important to test theoretical approaches experimentally in various ways. 

To overcome the aforementioned difficulties in the determination of \V u t,\, we proposed 
fi"3| a method of extracting \ V u b\ from inclusive charmless semileptonic decays of B mesons. 
We proposed to measure the lepton pair spectrum, i.e., the decay distribution of the kine- 
matic variable £ u = (q° + |q|)/Me i n the i?-meson rest frame, where q is the momentum 
transfer to the lepton pair and Mb denotes the B meson mass. Because most of B — > X u lv 
(£ = e or fi) events have a value of £ u beyond the threshold allowed for B — ► X c tv decay, 
£ u > 1 — Md/Mb = 0.65 with Mr> being the D meson mass, this kinematic requirement 
provides a powerful tool for background suppression. If nothing else, the decay distribution 
of £ u is of direct interest from an experimental standpoint. The sum rule can then be used to 
extract \V u h\ from the weighted integral of the measured £ u spectrum. The sum rule follows 
from the light-cone expansion for inclusive charmless semileptonic decays of B mesons and 
6-flavored quantum number conservation. The sum rule is thus independent of phenomeno- 
logical models. Moreover, the sum rule does not receive any perturbative QCD correction. 
Therefore, this method is not only experimentally very efficient, but theoretically also quite 
clean, allowing a precise determination of \V u b\ with a minimal overall (experimental and 
theoretical) error. 

The sum rule requires measuring the ^ u spectrum in inclusive charmless semileptonic B 
meson decays, weighted with £~ 5 . The purpose of the present paper is to refine the analysis 
of fpL3| . We calculate the perturbative QCD correction to the £ u spectrum to order a s . We 
investigate in detail how gluon radiation and hadronic bound state effects affect the shape 
of the £ u spectrum. 

The rest of the paper is organized as follows. In Sec. II the sum rule for inclusive charmless 
semileptonic decays of B mesons is briefly described. We discuss how to extract \V u b\ from 
experiment exploiting the sum rule. In Sec. Ill we study the £ u spectrum, including both the 
leading perturbative and nonperturbative QCD corrections. The exploration is extended to 
the weighted £ M spectrum in Sec. IV. Finally in Sec. V we present our conclusions. 



II. SUM RULE 

The sum rule for the decay distribution of £ u in inclusive charmless semileptonic decays 
of B mesons has been derived [|TB| from a nonperturbative treatment of QCD based on 



light-cone expansion. The large B meson mass sets a hard energy scale in the reactions, so 
that the light-cone expansion is applicable [f|-§| to inclusive B decays that are dominated by 
light-cone singularities, just as deep inelastic scattering. For inclusive charmless semileptonic 
decays of B mesons, 96% of phase space has q 2 > 1 GeV 2 . This dominance of the high-g 2 
region in the whole phase space < q 2 < M B renders the light-cone expansion especially 
reasonable. The theoretical expression for the weighted integral of the distribution of £ M 

r 1 1 dT 

S = / di u - — {B^X u iv) (1) 
Jo S d£ u 



is very clean. In the leading-twist approximation of QCD, it was found [13] that 
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S = \Vrt\ 2 ^^(B\1rf»b\B)-^ , (2) 
I ub \ 1927r 3 7 2M? ' 



where P is the momentum of the decaying B meson. The masses of leptons, the u quark and 
the 7r meson are neglected. Because ^-flavored quantum number conservation under strong 
interactions implies 

{B\b^b\B) = 2P», (3) 
the leading-twist contribution of nonperturbative QCD to the observable S in Eq. (0) is 



precisely calculable from first principles. The result is the sum rule [13 



S = f 1 d^LJL iB ^xju) = \V ub \ 2 ^l-. (4) 
Jo &d£u 192tt 3 v ; 

The conserved current Ir/^b is not subjected to renormalization and the corresponding 
anomalous dimension vanishes thanks to the Ward-Takahashi identity. Hence we recognize 
that the sum rule (f|) does not receive any perturbative QCD correction. 

It might seem surprising that one is able to calculate S so precisely in QCD when the 
experimentalist measures hadrons, and it is well known that hadrons cannot be understood 
completely within perturbation theory. There is a physical way to make clear that we 
really do not need to know anything about how hadrons are formed in order to predict S 
in inclusive B decays. The observable S measures something fundamental - namely the 
6-flavored quantum number carried by the B meson, and is basically determined by the 
underlying global symmetry of QCD. It is therefore insensitive to hadronic bound state 
effects. 

A measurement of the observable quantity S in inclusive charmless semileptonic decays 
of B mesons will lead to a determination of \V u b\ by use of the sum rule @ which contains 
no additional parameter. This determination of \V u b\ is free of perturbative QCD uncertain- 
ties. Meanwhile, the dominant hadronic uncertainty is avoided. Higher-twist corrections of 
nonperturbative QCD to the sum rule (f|) are suppressed by a power of Aq CD /M^. Once 
the observable S is measured, then Eq. @ yields |T4&| with small theoretical uncertainty. 
Moreover, this method is not only theoretically quite clean, but experimentally also very 
efficient in the discrimination between b —>■ u signal and b — > c background, as we shall 
discuss. 



III. THE £ u SPECTRUM 

In this section we describe the theoretical prediction for the £ u spectrum. We explore 
step by step perturbative and nonperturbative QCD effects on the £ u spectrum. We take 
into account first the perturbative QCD correction to order a s and then the leading non- 
perturbative QCD correction. 
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FIG. 1. The £ u spectrum in the tree- level and virtual gluon exchange processes b 
6-quark mass m;, = 4.9 GeV is used. 



utv. The 



A. Leading perturbative QCD correction 



Ignoring QCD corrections, the tree-level £ u spectrum in the free quark decay b 
the 6-quark rest frame takes the form 



r di u 



(b -> utv) = 5 ^ 



Mr 



utv in 



(5) 



with 



r n = 



G 2 F ml\V ub \ 2 

1927T 3 



(6) 



where m b represents the 6-quark mass. The resulting spectrum is a discrete line at £ u = 
itl^/Mb- This is simply a consequence of kinematics that fixes £ u to the single value m^/Ms, 
no other values of £ u are kinematically allowed in b — > utv decays. This is also the case for 
b — > utv processes with virtual gluon emission. Therefore, the spectrum in free quark 
decays is still a discrete line at = m^/M^, shown in Fig. 1, even if virtual gluon emission 
occurs. However, the above kinematic relation no longer holds for free 6-quark decays with 
gluon bremsstrahlung, and hence the spectrum expands downward below the parton-level 
end point = m^/M^. 

Including both virtual gluon emission and gluon bremsstrahlung, the differential decay 
rate as a function of £, u in the 6-quark rest frame is calculated to order a s to be 



1 dT 

To d^ u 



(b 



utv) = 5 
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(7) 



where r = 1 — Mb^/wi varying in the range < r < 1, corresponding to the kinematic 
range < £ u < nib/ Mb at the parton level of quarks and gluons. The distribution [g(r)] + 
is defined to coincide with the function g(r) for all values of r greater than 0, and to have 
a singularity at r = such that the integral of this distribution with any smooth function 
t(r) gives 



dr [g{r)]+t{r) = / dr g(r)[t(r) - t(0)]. 



(8) 



The result for the perturbative £ M spectrum with the QCD radiative correction^ to order ot s 
is shown in Fig. 2. The perturbative ^ u spectrum is singular at its end point £ u = m h /M B 
due to infrared divergences. We observe that gluon bremsstrahlung generates a small tail 
below £ u = mi,/ Mb- 

Integrating Eq. (0) over £ u yields the total perturbative decay rate 



r(6 -> utv) = r 



1 - 



2a. 
3vr 
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25 
T 



(9) 



This agrees with the well-known result obtained in [H|. 



B. Leading nonperturbative QCD correction 

To calculate the real physical decay distribution in inclusive decays B — * X u £u, pertur- 
bative QCD alone is not sufficient. We must also account for hadronic bound state effects 
due to the confinement of the b quark inside the B meson. We now consider nonperturbative 
QCD effects on the ^ u spectrum. In the framework of the light-cone expansion, the leading 
nonperturbative QCD effect is incorporated in the 6-quark distribution function || 

/(0 = ^^e^(5|6(0)^exp[^£^,(z)]%)|i?)| y2=0 , (10) 



1 The sum rule (||) does not receive any perturbative QCD correction because it results from the 
conserved current. The numerical calculation by integrating Eq. (]?]) with the factor £~ 5 has not 
reached the absolute zero but yielded a tiny correction (about 0.5%). We have checked that this 
is only a numerical artifact due to the limited accuracy of numerical integration by computer, 
especially given the rapid oscillation of the integrand near the end point £ u = 0. 
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FIG. 2. The perturbative £ u spectrum with the QCD radiative correction to order a s , obtained 
from input values m& = 4.9 GeV and a s = 0.2. 



where V denotes path ordering. The distribution function /(£) has a simple physical inter- 
pretation: It is the probability of finding a fe-quark with momentum £P inside the B meson 
with momentum P. The real physical spectrum is then obtained from a convolution of the 
hard perturbative spectrum with the soft nonperturbative distribution function: 



dT 



r 1 dT 



u£v,p b = £P), 



(11) 



where the 6-quark momentum pj, in the perturbative spectrum is replaced by £P. The an- 
alytic result for the perturbative spectrum to order a s is given in Eq. (0). The interplay 
between nonperturbative and perturbative QCD effects has been accounted for since con- 
finement implies that free quarks are not asymptotic states of the theory and the separation 
of perturbative and nonperturbative effects cannot be done in a clear-cut way. 

Equation (pd|) demonstrates that the physical £ M spectrum depends on the distribution 



function /(£). The definition of it, Eq. (ID), involving the 5-meson matrix element of the 
non-local 6-quark operators separated along the light cone, makes clear that the distribution 
function is a nonperturbative quantity. Although a complete calculation of the distribution 
function in QCD is impossible at present due to our ignorance of nonperturbative QCD, 
some basic properties of it are known ||. The distribution function is universal in the 
sense that the same distribution function also summarizes the leading nonperturbative QCD 
contribution in inclusive radiative B decays B — > X s 7. It is gauge invariant and obeys 
positivity. It has a support between and 1 and is exactly normalized to unity because 



7 



of 6-flavored quantum number conservation^. It contains the free quark decay as a limiting 
case with /(£) = 5(£ — m b jMB). In the free quark limit, Eq. (11) consistently reproduces 
the free quark spectrum. 

In addition, the mean /i and the variance a 2 of the distribution function were deduced 
H using operator product expansion and heavy quark effective theory (HQET) B: 



Mr 



1 + 



5E h 



(12) 



a 



( m b 



B 



2K h 



5E, ^ 2 



(13) 



where E b = K b + G b and K b and G b are the dimensionless HQET parameters of order 
(Aqcd/^6) 2 , which are often referred to by the alternate names Ai = — 2m 2 fT b and A 2 = 
— 2m 2 G b /3. The parameter \ 2 can be extracted from the B* — B mass splitting: A 2 = 
(M B * — M|)/4 = 0.12 GeV 2 . The parameter Ai suffers from large uncertainty. 

The mean value and variance of the distribution function characterize the location of the 
"center of mass" of the distribution function and the square of its width, respectively. They 
specify the primary shape of the distribution function. From Eqs. (|P2|) and (|13"D we know 
that the distribution function is sharply peaked around £ = // ~ m^/M^ close to 1 and its 
width of order Aqcd/Mr is narrow, suggesting that the distribution function is close to the 
delta function form in the free quark limit. 

Nonperturbative QCD methods such as lattice simulation and QCD sum rules could help 
determine further the functional form of the distribution function. The distribution function 
can also be extracted directly from experiments of inclusive semileptonic O] or radiative || 
decays of B mesons. The universality of the distribution function implies great predictive 
power: Once the distribution function is measured from one process, it can be used to make 
predictions in all other processes in a model-independent manner. Since these are as yet not 
done, we perform the calculations using the parametrization || of the distribution function 



/(0 = N 



£(i-0 c 



[(£ " a) 2 + b 2 } 



(14) 



where a, (3, a, and b are four parameters and iV is the normalization constant. The 



parametrization (|i4D respects all the known properties of the distribution function, in par- 
ticular the strong constraints of the sum rules (^) and (|T3"D. 

The £ u spectrum can be calculated using Eqs. ([TT|), ([^D and ([141). Including both the lead- 
ing nonperturbative and perturbative QCD corrections, we obtain the £ u spectrum shown in 
Fig. 3, using a s = 0.2, a = (3 = 1, a = 0.9548 and b = 0.005444. Here both values of a and 
P for the four-parameter distribution function (|TJ]) are preset to be 1, but in general they 
need not be integers. The values of a and b are then inferred from the sum rules (|12D and 



2 The same conservation law results in the sum rule 
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FIG. 3. The £ n spectrum including both the leading nonperturbative and perturbative QCD 
corrections. 



(|T3|) using nib = 4.9 GeV and Ai = —0.5 GeV 2 , giving the mean /i 
1 — 0.006 for the distribution function. 



0.93 and the variance 



a 



Bound-state effects lead to the extension of phase space from the parton level to the 
hadron level, also stretch the spectrum downward below m&/ Mb, and are solely responsible 
for populating the spectrum upward in the gap between the parton-level end point £ u = 
nib/ Mb and the hadron- level end point £ u = 1. The interplay between nonperturbative and 
perturbative QCD effects eliminates the singularity at the end point of the perturbative 
spectrum, so that the physical spectrum shows a smooth behavior over the entire range of 
£u, < £ u < 1, as m Fig- 3. Therefore, nonperturbative QCD effects play a crucial role in 
shaping the £ u spectrum. 

By integrating Eq. (|ll|) over £ u , we obtain the fraction of B — > X u lv events above the 
charm threshold allowed for the predominant B — > X c iv decays, defined as 



R=——^——f 1 d£ u ^(B 
T(B -> X u tv) Ji-m d /m b d£ u 



XJv) 



(15) 



We find that R = 79% for the spectrum shown in Fig. 3. This result refines our previous 
result fT3|. The earlier calculation without QCD radiative corrections found R = 99%. 



Gluon bremsstrahlung streches the £ M spectrum downward, and hence gives rise to a decrease 
of the fraction of B — > X u lv events above the charm threshold. Nevertheless, most of 
B — > X u iv events remain above the charm threshold. Thus the kinematic cut on the 
observable quantity £ u is very efficient in disentangling B — * X u iv signal from B — > X c lv 
background. This efficiency can be explained by the uniqueness of the £ u spectrum: The 
£ u spectrum stemming from the quark-level and virtual gluon exchange processes would 
only concentrate at £ u = m b /M B , shown in Fig. 1, solely on kinematic grounds, and gluon 
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bremsstrahlung and hadronic bound state effects smear the spectrum about this point, but 
most of the decay rate remains at large values of £ u , as revealed by Figs. 2 and 3. 

The fraction R of course depends on forms of the distribution function. However, we 
find that R is relatively insensitive to forms of the distribution function once the mean and 
variance of it, which are known from HQET, given by Eqs. flT2| ) and (|T3|), are kept fixed. 
Therefore, the above calculation of R can be considered as a typical estimate of the fraction 
of B — ► Xjjiu events above the charm threshold. 



IV. THE WEIGHTED SPECTRUM 

The weighted spectrum ^~ 5 dT(B —>■ X u £v) / d£ u is more directly relevant for the mea- 
surement of the observable S in inclusive charmless semileptonic decays of B mesons. Mea- 
surements of S can be obtained from an extrapolation of the weighted spectrum measured 
above the charm threshold to the full phase space available in B — > X u iv decays. While the 
normalization of the weighted spectrum given by the sum rule @ does not depend on the 
6-quark distribution function /(£), thus being model-independent, the shape of the weighted 
spectrum does. In this section we perform a detailed analysis of the shape of the weighted 
spectrum, using Eqs. flTT|), (|7]) and (Hp. 



We first explore the impact of gluon radiation on the shape of the weighted £ u spec- 
trum. We show in Fig. 4 the weighted £ u spectram without and with the QCD radiative 
correction. It is evident from Fig. 4 that the shape of the weighted £ u spectrum receives a 
significant correction due to gluon radiation. However, the shape of the spectrum appears 
to be insensitive to the value of the strong coupling a s , varied within a reasonable range. 

We investigate next the sensitivity of the spectrum to the form of the distribution func- 
tion, which would reflect the impact of hadronic bound state effects. For this purpose we 



choose to use two very different forms [15| for the distribution function, albeit having the 
same mean value and variance. The calculated weighted spectra are shown in Fig. 5, taking 
into account the QCD radiative correction to order a s . The weighted spectrum exhibits a 
strong dependence on the form of the distribution function, even though the mean value 
and variance of the distribution functions are the same. 

Generally, since the quark-level processes, exclusive of gluon bremsstrahlung, generate a 
discrete line, the shape of the £ u spectrum directly reflects the inner long-distance dynamics 
of the reaction. This argument elucidates the strong variation of the weighted spectrum 
with the form of the distribution function, illustrated in Fig. 5. 



In fact, ignoring QCD radiative corrections, one obtains from Eqs. (11) and (|5]) the 
weighted spectrum in the B rest frame 

It makes clear that the resulting weighted spectrum is directly proportional to the distri- 
bution function. In other words, the weighted £ u spectrum is most sensitive to the distri- 
bution function. This salient feature renders the weighted spectrum idealy suited for the 
direct extraction of the distribution function from experiment. After an observed spectrum 
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FIG. 4. The weighted £ u spectrum without (dotted curve) and with (solid curve: a s = 0.2, 
dashed curve: a s = 0.27) the QCD radiative correction. The four parameters for the distribution 
function are taken to be a = (3 = 1, a = 0.9548, and b = 0.005444. All the spectra are normalized 
to have unit area. 

^~ 5 dT(B — > X u iv)jd^ u is radiatively corrected, it is nothing but the 6-quark distribution 
function if higher-twist terms are neglected. 

Extrapolation of the weighted £ u spectrum measured above the charm threshold to the 
full £ u range requires a theoretical calculation of the rate ratio 



1 r 1 
W=- 

b Jl-M D /M E 



1 dT 

WnWu 



B -> XM 



(17) 



where S is denned in Eq. (|l|). We calculate the ratio W using Eqs. (§), flTl"D , (0) and 
and the theoretical uncertainties are estimated as follows: 

We study the variation of W with respect to the mean value and the variance of the 
distribution function setting a — [3 — 1 in Eq. (|H|). Actually, this amounts to the study of 
the ratio W as functions of and Ai, since, essentially, the mean value of the distribution 
function is determined by the 6-quark mass and its variance is determined by Ai according 
to the sum rules in Eqs. ( |I2D and (|T3|). At present, the estimated values of the 6-quark mass 
and Ai vary in the ranges 



m b = 4.9 ±0.15 GeV, 
Ai = -(0.5 ±0.2) GeV 2 . 



(18) 
(19) 



The variation of m& leads to an uncertainty of 8% in W if other parameters are kept fixed. 
A small uncertainty of 3% in W results from the variation of Ai- In other words, the ratio W 
displays a strong dependence on the mean value of the distribution function of the b quark 
inside the B meson, but is insensitive to the variance of the distribution function. 
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FIG. 5. Comparison of the weighted £ u spectra in two forms of the distribution function. 
The solid curve corresponds to the four-parameter distribution function (|l4|) with a = (3 = 1, 
a = 0.9548, and b = 0.005444, while the dashed curve corresponds to the four-parameter distri- 
bution function ( |l~4|) with a = (3 = 2, a = 0.9864, and b = 0.02557. These two forms of the 
distribution function have the same mean value fi = 0.93 and variance a 2 = 0.006. The strong 
coupling is taken to be a s = 0.2. Both spectra are normalized to have unit area. 

We examine the further sensitivity of the rate ratio W to the form of the distribution 
function when keeping the mean value and variance of it fixed, by varying the values of 
the two additional parameters a and f3 in the parametrization ([14]). We estimate that the 
variation of W is 3% if the form of the distribution function is changed but with the same 
mean value and variance. 

We estimate the uncertainty due to the truncation of the perturbative series in Eq. (0) 
by varying the renormalization scale between rrib/2 and 2m&. We find that an uncertainty 
of 8% in the ratio W stems from the renormalization scale dependence. 

This analysis implies that at present the theoretical error in the calculation of W has 
two main sources: the value of m& (or equivalently, the mean value of the distribution 
function) and the renormalization scale dependence. Finally, adding all the uncertainties in 
quadrature we arrive at 

W = (76 ± 9)% (20) 

with an error of 12%. 

V. CONCLUSIONS 

We have studied the £ u spectrum in inclusive charmless semileptonic decays of B mesons. 
The perturbative QCD correction to the spectrum is calculated to order a s . The leading 
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nonperturbative QCD effect is calculated using light-cone expansion and heavy quark effec- 
tive theory. The £ u spectrum is unique in that the tree-level and virtual gluon exchange 
processes b — > utv at the parton level generate a trivial £ u spectrum - a discrete line at 
£ n = m^/Ms ~ 0.93 in the 6-quark rest frame, which is well above the charm threshold 
£ n = 1 — Md/ 'Mb = 0.65. Gluon bremsstrahlung results in a small tail in the lepton pair 
spectrum below mj/ Mb- Bound-state effects lead to the extension of phase space from the 
parton level to the hadron level, also stretch the spectrum downward below nib/ Mb, and are 
solely responsible for populating the spectrum upward in the gap between the parton-level 
endpoint £ u = m h /M B and the hadron- level endpoint £ u = 1, smoothing out the singularity 
in the perturbative spectrum. As a result of these two distinct effects, the lepton pair spec- 
trum spreads over the entire physical range < £ u < 1. Still, we find that about 80% of 
the lepton pair spectrum in B — > X u lv lies above the charm threshold, £ u > 1 — Mb /Mb- 
This kinematic cut is most efficient in the suppression of the background of B semileptonic 
decays into charmed particles. 

The observable S can be measured by the extrapolation of the weighted spectrum 
£,~ 5 dT(B — > X u tv)jd^ u measured above the charm threshold to the entire phase space. 
Gluon bremsstrahlung and hadronic bound state effects strongly affect the shape of the 
weighted £ u spectrum. However, the shape of the weighted £ u spectrum is insensitive to the 
value of the strong coupling a s , varied in a reasonable range. The overall picture appears 
to be that the weighted £ u spectrum is peaked towards larger values of £ u with a narrow 
width. The contribution below £ u = 0.65 is moderate and relatively insensitive to forms of 
the distribution function. This suggests that extrapolating the weighted £ u spectrum down 
to low £ u would not introduce a considerable uncertainty in the value of 5". Quantitatively, 
our analysis determines the rate ratio for the extrapolation of the weighted spectrum to 
be W = (76 ± 9)% with an error of 12% at present. 

Another interesting use of the weighted £ u spectrum is its utility for directly extracting 
the 6-quark distribution function. The universality of the distribution function implies that 
the distribution function extracted from inclusive semileptonic B decays can be used to 
make predictions in inclusive radiative decays B — > A s 7 in a model-independent manner 
and vice versa. In particular, a measurement of the distribution function from the photon 
energy spectrum in B — ► X s ^j || would be very useful to improve the measurement of S 
from extrapolation of the weighted £ u spectrum. 

An improved knowledge of the form of the distribution function is important for an error 
reduction in extrapolation. More precisely, if by direct measurements or theoretical studies 
the uncertainty in the mean value of the distribution function can be improved from currently 
3% to 1.5% in the future, the resulting uncertainty in the ratio W for the extrapolation would 
decrease from 8% to 3%. A calculation of the 0(a 2 s ) perturbative QCD correction to the £ u 
spectrum is also important for improving the accuracy of the extrapolation. 

The physical observable quantity S is connected with \V u b\ via the sum rule (f|) for 
semileptonic B decay in the leading-twist approximation of QCD. The sum rule is a funda- 
mental, model-independent prediction of QCD. Note that the sum rule does not rely on the 
heavy quark effective theory. No arbitrary parameter other than \V u b\ enters the sum rule. 
The dominant hadronic uncertainty is avoided in the sum rule. As important, there is no 
perturbative QCD modification of this sum rule, so that the potential source of theoretical 
uncertainty associated with perturbative QCD calculations is totally averted. 
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The kinematic cut on £ u , £ u > 1 — Md/Mb, and the semileptonic B decay sum rule, 
Eq. (|j), offer an outstanding opportunity for the precise determination of \V u b\ from the ob- 
servable S. We wish to emphasize that this method is both exceptionally clean theoretically 
and very efficient experimentally in background suppression. 

There remain two types of theoretical uncertainties in the determination of \V u b\. First, 
higher-twist (or power suppressed) corrections to the sum rule cause a theoretical error of 
order Aq CB /M b ~ 1% on \V u b\. A quantitative study of higher- twist effects could further 
reduce this small theoretical uncertainty associated with simply extracting the value of \V u b\ 
from the measured value of S. Second, the extrapolation of the weighted £ u spectrum to low 
£ u gives rise to a systematic error for the measurement of S. The status of this additional 
theoretical uncertainty associated with extrapolation, as well as how to improve it, have 
been discussed above. 

Eventually, the error on \V u b\ determined by this method would mainly depend on how 
well the observable S can be measured. To measure S experimentally one needs to be able to 
reconstruct the neutrino. This poses a challenge to experiment. Given the unique potential 
of determining \V u b\, we would urge our experimental colleagues to examine the feasibility 
of the method. 
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